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The bound state in the system of three protons and an electron (pppe) under a homoge- 
neous strong magnetic field where the protons are situated in the vertices of an equilateral tri- 
angle perpendicular to the magnetic field lines is found. It is shown that for magnetic fields 
B = 10^^ - 4.414 X lO^^G the potential energy curves as a function of the internuclear distance 
R have an explicit minimum. For all magnetic fields studied, the binding energy of the triangular 
configuration is less than the binding energy of the linear parallel configuration (A. Turbiner et al. 
JETP Lett. 69, p. 844). In the contrary to the linear case, the binding energy decreases with a 
- , magnetic field growth, while the equilibrium internuclear distance slowly increases. 

, Recently a quantitative study of the system of three protons and an electron {pppe) in a strong magnetic field gave 
■ theoretical evidence of the existence of a bound state for magnetic fields B > lO^^G |jl[. It manifests the existence of 
(— ( the exotic molecular ion H^'^^\ Many years ago it was indicated by Kadomtsev and Kudryavtsev [||, and Ruderman 
d the possible existence of linear exotic molecular systems in presence of a strong magnetic field. Due to the enormous 

, Lorentz force acting on the electronic cloud, it shrinks in the direction transversal to the magnetic field lines leading to 
^\ ' an effective quasi-one-dimensionality of the studied systems. This makes the electron-nuclei attraction more effective 
T-H to compensate the Coulombic repulsion of nuclei. Such a fact suggests the possibility of the existence of exotic systems 
which do not exist without a magnetic field. It has been shown some time ago that a strong magnetic field can lead 
to the formation of linear hydrogenic chains iJ^, n > 2 situated along magnetic lines More recently the first 

quantitative study about the possible existence of an exotic bound state in the system of three protons and an electron 
{pppe) in a strong magnetic field was carried out 0| and it confirmed existence of H^^'^ ion. 
^nJ Although the linear parallel configuration (where protons are situated on a line parallel to the magnetic field line) 

for the system {pppe) seems to be optimal for very strong magnetic fields, other possible configurations have to be 
, considered in order to have a better understanding of the properties of such systems. 

' In this article we explore the possibility of the existence of the molecular ion H^^^ , where the protons are situated 
in the vertices of an equilateral triangle perpendicular to a homogeneous magnetic field. We limit our present 
investigation to magnetic fields B = 10^^ — 4.414 x lO^^G for which a non-relativistic approach is still valid (for 
a detailed discussion see [Q and references therein). Our study is restricted to the ground state within the Born- 
O Oppenheimer approximation. 

^ The present calculations are carried out in the framework of a variational method. Trial functions are chosen 

^ in accordance with the following criterion (see j|] — [0): (i) the trial function '^t{x) should include all symmetry 
. . properties of the problem in hand; (ii) the trial function for the ground state should not vanish inside the domain 

^ ' where the problem is defined; (iii) the potential Vt{x) = for which the trial function is an exact eigenfunction, 

,; should reproduce the original potential behavior near singularities as well as its asymptotic behavior. This prescription 
^ ■ has been successfully applied to the study of the H2 ion in strong magnetic fields giving the more accurate results for 
5^ B > S, as well as to the study of the linear parallel configuration of the h'^^^ and Hf'^'^ ions in a strong magnetic 
fields 

We consider a system of one electron and three identical infinitely-heavy centers of unit charge situated on the 
vertices of an equilateral triangle with internuclear distance R (which is the length of the triangle side) lying in a 
plane perpendicular to the magnetic field of strength B directed along the z axis, B = (0, 0, B). 

The potential corresponding to the system we study is given by 
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where the quantity p = ya^+y^ is the distance from the electron position to the z-axis, ^1^2, 3 are the distances 
from the electron to the first, second and third centers, respectively, R is the distance between centers (see Fig. |l| 
for notations) and is the z-component of the angular momentum operator. The linear Zeeman effect term in the 
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potential {B£z) can be dropped off since we use real trial functions Q (see below). Spin Zeeman effects are neglected. 
Through the paper the Rydberg is used as the energy unit. For the other quantities standard atomic units are used. 

The trial functions we are going to use for the present problem are similar to those were exploited in for the 
study of the linear parallel configuration. These functions were built according to the criterion described in ||^ — 
. They contain the basic features of Coulomb systems in a magnetic field, as well as the symmetry property under 
permutation of the three charged centers. In particular, it implies that the potentials corresponding to those functions 
reproduce the Coulomb- like behavior near the centers and the two-dimensional oscillator behavior in the (x, y) plane 
at large distances (see below). 

The simplest of those functions is a function of the Heitler-London type multiplied by the lowest Landau orbital: 



Vj/^ _ g-ai(ri+r2+r-3)-^i-Bp^/4^ j^2) 

(cf. Eq. (2) in Ref. [Q), where ai and /3i are variational parameters. Since we consider the distance between centers as 
an extra variational parameter, it has in total three variational parameters. This function gives an adequate description 
of the covalent coupling of the system near equilibrium. As illustration we show the potential corresponding to this 
function: 

Vx = = ia\ ~ l3iB - 2ai y - + 2al V • + ai/JiSV ^— ^ + ^ ' , 3 

i—1 i<j i—1 

■hi ■ hj = ^— [(x - Xi){x - Xj) + (y - yi)(y - + {z - Zi){z - Zj)] , 

where fii (i = 1,2,3) is the unit vector in the direction of the vector pointing from the position of the i-th center to 
the position of the electron. 

A second trial function is a Hund-Mulliken type function multiplied by the lowest Landau orbital: 



(cf. Eq. (4) in Ref. ||i|). Here a2, P2 are variational parameters. This function describes a ionic coupling of a hydrogen 
atom with two charged centers. 

Another suitable trial function which describes a ionic coupling between an H2 ion and a proton is given by 



\[>g — ^g-a3('-i+»-2) _|_ g-aaCri+'-a) _|_ g-a3('-2+'-3)~j g-/53Sp^/4 ^ j-g-j 

(cf. Eq. (5) in Ref. Q]) where a3,/33 are variational parameters. Since the H+ ion has the lowest total energy among 
the one-electron systems for B < lO^^G an important contribution coming from this trial function is expected. 

In order to include in a single trial function the different physical behavior near equilibrium and at large distances, 
appropriate interpolations of the trial functions (^, (^ and (^) are done: 

(i) A natural interpolation is given by a non-linear superposition of the form: 




(cf. Eq. (6) in Ref. §) where a4, Qfs, ae and are variational parameters, and the sum is over all permutations 
of the parameters {a^, as, ag}. 

If all parameters coincide (04 = (15 =05 — cei): the function (^ reduces to the Hcitlcr-London type function 
(0) . Function m) also reduces to the Hund-Mulliken type wave function (H) when only one parameter is non-zero 



^For a real trial function the expectation value of the z-component of the angular momentum {l^} equals to zero and leads to 
a vanishing contribution to the variational energy emerging from the linear Zeeman effect term . 
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say a4 = and 05 = ag = 0. If two parameters are non-zero, and equal, say, = =03, and ae = it 
reduces to the trial function (^. When all parameters are different among them and different from zero the 
function provides us with a 3-center modification of a Guillemin-Zener type function. The function ^ has 
in total five variational parameters (including the internuclear distance R as one of the parameters). 

(ii) Another more immediate interpolation is given by a linear superposition of the functions (0), (H) and (p|). 



*5_i, = Ai1'i+^2*2+^3*3, (7) 

where Ai, A2, are taken as extra variational parameters. Since ^1,2,3 are not orthogonal, the parameters 
^1,2.3 do not have the usual meaning of weight factors. 

For the present calculations we use a linear superposition of the above interpolations (^), given by: 

^-6 = Ai^-i + ^2^-2 + ^3*3 + A4^'4-„is, (8) 

where Ai, A2, A3, A4 are considered again as extra variational parameters with no meaning of weight factors. The 
function combines in a single function form, a suitable description of the system for the different physical regimes 
(near equilibrium and large internuclear distances for strong and very strong magnetic fields). Therefore we guess it 
should provide a relevant approximation for the ground state of the system. This was indeed the case of the linear 
parallel configuration of H^^^ (see Ref. ||l[). This function has in total fourteen variational parameters. 

The variational procedure is carried out using the standard minimization package MINUIT from CERN-LIB. We 
use the integration routine DOIFCF from NAG-LIB. All integrals are calculated with relative accuracy ^ 10~^. 

The results of our calculations are presented in Table I. For all considered magnetic fields B = 10^^ — 4.414 x lO^'^G 
we found a minimum in the potential energy curve as a function of the internuclear distance R (see below Fig. 3) 
indicating the formation of a bound state. For such magnetic fields the total energy of the triangular configuration 
is always larger than the total energy of the corresponding linear parallel configuration. Indeed this shows that 
for very strong magnetic fields the linear configuration is the most favored. However, a new striking unexpected 
phenomenon appears when one considers a triangular configuration. While for the linear configuration the binding 
energy^ increases and the internuclear distance decreases as the magnetic field grows, the opposite is true for the 
triangular configuration. Namely, the binding energy is slowly decreasing and the internuclear distance is slowly 
increasing (!) both being almost constant. A similar behavior occurs with the longitudinal localization length of the 
electron (|z|) which also slowly increases in the triangular configuration with a magnetic field growth, opposite to 
what happens in the linear configuration where it decreases rather sharply. The transversal size of the electronic cloud 
{p) is very close to the corresponding cyclotron radius i3~^/^(a.u.) for both cases. We do not have a clear physical 
picture of this phenomenon and it needs further consideration. 

A surprising feature is revealed if we study different trial functions. In Table II it is shown the variational results for 
the total energy Et for two quite arbitrarily chosen values of the internuclear distance, R = 0.45 and R = 2.527. For 
magnetic fields 10"'^^ ^ B ^ 10^^ G an adequate description of the system for internuclear distances R ~ 0.45 is given 
by the trial function (|^) modelling a ionic coupling, -t-proton. In this case, the single function (^ gives the lowest 
total energy among the trial functions 'I'i,2,3,4- On the other hand, the trial function (g) modelling a covalent coupling 
gives the lowest total energy for R ~ 2.527 giving a hint that a covalent coupling dominates when the internuclear 
distance is large compared to the equilibrium distance. Amazingly, for B ^ lO^^G the variational energy seems to be 
insensitive to the specific form of the trial function (!). 

The electronic density distribution jvPl^ in the plane z = is shown in Fig. 2. For all magnetic fields B — 
10^^ — 4.414 X 10"'^^ G the electronic density distribution j^fp exhibits a single maximum located at the center of 
the triangle {x — 0,y — 0). As the magnetic field grows the electronic cloud becomes more concentrated around 
X = 0, y = 0, while its distance to the protons is slowly increased. 

(2+) 

From our results we can draw the conclusion that the exotic molecular ion can exist in the presence of a 

strong magnetic field in a triangular configuration of the charged centers which is unstable towards a decay mode to 



■^The binding energy is defined as the affinity to keep the electron bound Et = B — Et- B is given in Ry and thus has a 
meaning of the energy of free electron in magnetic field. 
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the linear configuration. In particular, in the present article we have shown that the molecular ion i/g can also 
exist in the equilateral (fixed nuclei) triangular configuration perpendicular to the magnetic field. 

The results of our variational study for the ground state show that, contrary to the case of the linear parallel 
configuration, the binding energy of the system decreases, and the internuclear distance increases, as the magnetic 
field grows from B — lO^^G up to 4.414 x lO^^G. A more extended study of this as well as other configurations is 

(2+) 

needed in order to have a better understanding of the properties of the exotic system iJg . In particular, a proper 
contribution coming from the linear Zeeman effect term in (m) has to be taken into account. However, the fact that the 
linear Zeeman effect contribution vanishes in our consideration (which occurs when real trial functions are considered) 
looks relevant in the region of very strong magnetic fields B = 10^^ — 4.414 x lO^'^ G where this term can be neglected 
anyhow. 

Although in general it is not clear so far (see discussion in 0) how an adiabatic separation of the electronic and 
the nuclear motion can be performed in the presence of a strong external magnetic field, our fixed-nuclei approach is 

(2+) 

a good starting point to study of the exotic molecular ion iJg 
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TABLE I. Variational results for the triangular configuration of the ion in a strong magnetic field. The total energy 

Et and the binding energy Eh = B — Et, are in Ry. The equilibrium internuclear distance Req and the expectation values 
of the longitudinal and transversal localization length of the electron {\z\) and (|pj) are in a.u.. The corresponding quantities 
for the fff ion in the linear parallel configuration were taken from The energy of free electron and its cyclotron radius 
B~^^^ are also shown. 



B (Gauss) 


Et 


Eb 




<\z\> 


<P> 


Triangular Configuration 




40.5591 


1.985 


0.6237 


0.561 


0.192 




423.722 


1.719 


0.6414 


0.595 


0.061 




4252.724 


1.690 


0.6645 


0.610 


0.019 


4.414 X 10" 


18777.31 


1.674 


0.6757 


0.622 


0.009 


Linear Configuration 


10^^ 


36.429 


6.1151 


0.803 


0.864 


0.186 


10" 


410.296 


15.144 


0.346 


0.438 


0.060 


10" 


4220.09 


34.324 


0.165 


0.242 


0.019 


4.414 X 10" 


18723.88 


55.103 


0.110 


0.168 


0.009 


Free electron 


10^^ 


42.5441 








0.153 


10" 


425.441 








0.048 


10" 


4254.41 








0.015 


4.414 X 10" 


18778.98 








0.007 



TABLE II. Comparison of the total energy Et for the triangular configuration of the system (pppe) for two arbitrary 
internuclear distances R — 0.45 and R — 2.527 obtained with trial functions Vl/1,2.3,4, and with their linear superposition given 
by trial function 'I'e. 



B (Gauss) 


R 




*2 


*3 


*4 


*6 


10" 


0.45 


41.10121 


41.12298 


41.09077 


41.09591 


41.08986 


2.527 


41.58129 


41.58405 


41.58146 


41.58131 


41.58122 


10" 


0.45 
2.527 


423.9111 
424.4875 


423.9195 
424.4877 


423.9115 
424.4876 


423.9111 
424.4876 


423.9110 
424.4865 


10" 


0.45 


4252.956 


4252.956 


4252.956 


4252.956 


4252.956 


2.527 


4253.461 


4253.461 


4253.461 


4253.461 


4253.460 


4.414 X 10" 


0.45 
2.527 


18777.53 
18778.03 


18777.53 
18778.03 


18777.53 
18778.03 


18777.53 
18778.03 


18777.53 
18778.03 
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FIG. 1. Triangular configuration for the (pppe) system in a magnetic field B directed along the z axis. The protons are 
situated in the {x,y)-p\a,ne and their positions are marked by black circles. 
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FIG. 3. Potential energy curves as function of the internuclear distance R for the H.;^^^ molecular ion in the triangular 
configuration for magnetic field: (a) B = 10"G, (b) B = lO^^G, (c) B = The top of the energy scale corresponds to 

the value of the energy of free electron Ee = B. 



42.54 




423.5 



4254.4 - 

4254 
E(Ry) 
4253.5 

4253 - 

4252.5 
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